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Nomenclature 


A area 

AN annular mixer abbreviation 

AN-Stag annular mixer geometry with staggered orifice configuration inner/outer wall 

A,,, duct cross-sectional area, also A tot , in 2 

ACd also Aj, effective orifice area, in 2 

Ar area ratio (jet/mainstream) = = ACd/A m 

AMIX area weighted overall T or (() deviation from stoichiometry, Eq 6 

B area determined half width of the distribution function 

CIR circular mixer abbreviation 

d diameter of the orifice, inch 

Do outer diameter of the mixer, inch 

DR density ratio (jet/mainstream) 

DP/P total pressure loss across the mixing wall, % 

EINOx global NOx El 

f non-dimensional equivalence ratio, Eq 5 

f/a fuel to air ratio 

HO configuration nomenclature (hole-optimization) 

j radial vector direction 

J momentum-flux ratio (jet/mainstream) = M 2 /DR, also (MR) 2 /[(DR)(Aj/An,) 2 ] 

k tangential vector direction 

L orifice length, inch (longest dimension of the slot) 

m tot overall mixer mass flowrate, lbm/sec = m jet + m mam 

M mass flux ratio (jet/mainstream) = DR Vj e t AJ ma in 

MMIX mass flow weighted overall T or (J) deviation from equilibrium, Eq 7 

MR mass flowrate ratio (jet/mainstream) 

P total pressure, atm 

Pen r/R non-dimensional jet penetration from the outer wall of the mixer 
R radius of the mixing section, inch 

r radial distance from the centerline of the mixer, inch 

S orifice spacing in the circumferential direction, inch 

T temperature, R 

T.E. x/R orifice trailing edge axial location (for the circular mixer) 

T.E. 2x/H orifice trailing edge axial location (for the annular mixer) 

u local axial velocity, ft/sec 

U ma in approach mainstream axial velocity, ft/sec 

VR velocity ratio (jet/mainstream) = V je« / U ma j n 

V jet radial velocity of the jet, ft/sec 

W orifice width, inch, (shortest dimension of the slot) 

x axial distance from the leading edge of the orifice, inch 

p fluid density, lbm/ft 3 

<|> equivalence ratio (f/a) local / (f/a) stoi 
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Subscripts 


equil 

equilibrium 

j 

jet 

LZ 

lean-zone 

m 

mainstream, also (main) 

RZ 

rich-zone 

stoi 

stoichiometric 
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Introduction 


The design and development of gas turbine engines for the aeropropulsion and ground based 
power generation systems has been directed towards decreasing the gaseous emissions without 
adversely affecting the system performance. The environmental effects of carbon monoxide 
(CO), oxides of nitrogen (NOx), and unbumed hydrocarbon fuels have been investigated 
extensively for many years and their potential hazards are increasingly becoming a more world 
wide sensitive issue. In specific, the effect of the NOx emissions on the environmental structure 
of our planet demands close attention due to the present ozone layer depletion discoveries. 

The advancement in the gas turbine engine design has been channeled towards increasing 
both engine pressure ratio and rotor inlet temperature levels in order to increase the overall 
thermodynamic cycle efficiency thus reducing the specific fuel consumption. However, this 
effort has been adversely affecting both engine durability and gaseous emissions, especially that 
of NOx which is generally directly proportional to the cycle temperature in air breathing engines. 
The production of oxides of Nitrogen (generally classified into both NO and NO 2 gas emissions) 
are generally a function of the combustion system design, hot section temperature distribution, 
residence time, and localized fuel/air mixture. The NOx emission in a gas turbine engine can be 
divided into two main types: a) prompt NOx, which is generally associated with low tempera- 
ture and fuel-rich mixtures and generally occurring early in the combustion typically near the 
fuel injector and primary zone of the combustion liner, and b) thermal NOx which is dominated 
with the available dissociated Oxygen atoms from the air mixed with the fuel. The thermal NOx 
is generally the highest NOx production in a typical gas turbine engine specifically in non- 
premixed diffusion flame combustion systems. 

The development of an efficient and low-emission high speed civil transport combustion 
system demands increased insight in combustion chemical kinetics, efficient air/fuel mixture, and 
advanced materials and cooling techniques. The combustion liner cooling air flowrate 
significantly affects the combustion process by quenching the flame at localized regions near the 
wall thereby producing CO and also producing localized hot spots which increase the production 
of NOx. Current gas turbine combustion liner technology employs a single-stage combustion 
process in which both fuel and air are admitted into a controlled mixing chamber and are allowed 
to react. With the advent of higher temperature operating conditions, the reduction of NOx 
becomes a very difficult task to accomplish using single axial staged combustion. Therefore, 
alternative combustion methods are being explored by the main gas turbine engine manufacturers 
and other research organizations. Two of the main low NOx designs being developed are axially 
multistage combustion rich-burn/quick-mix/lean-burn (RQL), and lean premixed prevaporized 
(LPP). Both have advantages and disadvantages regarding operation range and hardware 
complexity. 

In this report, the advantages and limitations of the RQL design are being explored by 
focusing on the performance of the quick-mix region where generally most of the thermal NOx 
is produced. The successful performance of the RQL combustor depends on a quick and 
efficient mixing of the rich zone combustion products with compressor discharge air to 
effectively reduce the overall and localized equivalence ratio from about 1.8 to about 0.5. This 
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process must be accomplished with a minimum transient time (at or near equivalence ratio of 1) 
where most of the NOx is produced due to the high resulting temperature levels and Oxygen 

availability. 

An experimental and analytical effort (e.g. Bain, Smith, and Holdeman, 1992, .1993, and 
1994; Doerr and Hennecke, 1993; Doerr, 1994; Doerr, Blomeyer, and Hermecke, 1995, Hatch, et 
al. 1992a and 1992b; Hatch, Sowa, and Samuelsen, 1994; Howe, et al., 1991; Kroll, etal., 1993, 
Kroll, Sowa, and Samuelsen, 1995; Liscinsky et al., 1992; Liscinsky, Vranos, and Lohmann, 
1993; Liscinsky, True, and Holdeman, 1993, 1994, and 1995; Oechsle Mongia and Holdeman 
1992 1993 and 1994- Oechsle, 1995; Oechsle and Holdeman, 1995; Smith, Talpallikar, and 
Holdeman, 1991; Sowa, et al., 1994); Talpallikar, et al., 1991; Vranos, et al., 1991; and Zhu and 
Lai, 1992) is underway to study and identify the critical design and flow parameters affecting the 
mixing effectiveness and its impact on NOx formation. 


The overall scope of this program is to analyze several mixing configurations numerically 
and assess the major design parameters in order to apply this technology for the design of future 
RQL combustion systems for use in a high speed civil transport. This program is a 5 year NAS 
funded program in which over 100 different mixing configurations have been numerically 
analyzed spanning the domain of both non-reacting and reacting flowfields and both circular and 
annular mixing configurations. 


In the first stage of the program, 21 different non-reacting cases were summarized using a 
circular mixing geometry. This study by Oechsle, Mongia, and Holdeman, (1992) indicates that 
mixing flowfield significantly depends of the operating conditions. Several orifice types and 
operating conditions were investigated in this study. The mixing uniformity of the flowfield also 
is significantly changed by the penetration of the jets. Large unmixed regions were found behind 
the orifice. This was a liability in a circular mixer design which tended to produce rich pockets 
of mixture and high quantities of NOx. The study in Oechsle, Mongia, and Holdeman, (1992) 
was performed in a purely non-reacting environment 

Following the 1992 study, a second study by Oechsle, Mongia, and Holdeman, (1993) was 
made to compare the predicted results with experimental data obtained by the University of 
California at Irvine by Hatch, et al. (1992a and b). This parallel numerical-experimental study 
allowed a detailed comparison of both methods to determine their strengths and weaknesses. 
The results show that the general trends are similar for the experimental and numerical data, 
however, mixedness was generally more developed in the experimental data and steeper 
gradients were more characteristic of the numerical results. These results also showed that the 
effect of orifice geometry on the mixing flowfield and jet penetration for the measured results 
compared favorably to the numerically obtained values. It was also established that the mixing 
flowfield in both the numerical and experimental results would be inherently different since the 
initial conditions for both are different. In the numerical simulation, the jet and mainstream mass 
injections are uniformly distributed throughout the entire injection area. Velocity, temperature, 
and mass fraction are uniform across the entire injection area. This was not obtained physically 
in the experimental setup. It is also important to consider that the inlet turbulence levels also 
significantly affect the mixing profile, and though these levels were not measured in the 
experimental setup, it is assumed that they are generally different than the standard values used 
in the numerical approach. Furthermore, more accurate comparisons of numerical and 
experimental results necessitate the need to clearly quantify and reproduce as much as possible 
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the initial conditions obtained in the experimental setup before duplicating them in the numerical 
model since it is believed that the reverse will be more difficult to accomplish. 

In the following phase of the program, a comparison of the mixing flowfields in both 
reacting and non-reacting environments are shown in Oechsle, Mongia, and Holdeman, (1994). 
These results indicate good agreement in the conserved scalar fields in both reacting cases and 
non-reacting cases with similar orifice configurations and operating conditions. It was also 
evident from these results that the pockets of stoichiometric mixture near the outer wall behind 
the orifice were present in both reacting and non-reacting numerical calculations which verifies 
the results of the previous two studies. The statistical mixing non-uniformity parameters 
correlated with the observational data in predicting the best mixing configuration. 

In the following study Oechsle, and Holdeman, (1995), eighteen different circular mixing 
configurations were analyzed in a fully reacting environment at actual engine altitude cruise 
operating condition. The results indicate that the mixing uniformity, as defined statistically, did 
not correlate with NOx production trends. The expected results indicates that better mixing of 
the mainstream flow and crossflow jets would yield an increase in temperature uniformity which 
consequently leads to lower NOx production. These expectations are however skewed by the 
way the mixing uniformity is quantified. In the above mentioned study, the mixing non- 
uniformity depends on the deviation from the equilibrium equivalence ratio value. This 
definition does not consider the deviation from the stoichiometric equivalence ratio value where 
most of the NOx is presumed to be produced. The results also indicate that the round orifices 
produce less NOx at low J conditions (J=25) and that slanted slots produce less NOx at higher J 
conditions (J=52, 80). Therefore, jet penetration which is significantly affected by the orifice 
shape has a strong impact on the NOx production. 

In the present study, a 3-D numerical CFD code is used to predict the performance of an 
RQL mixing section. The mixing flowfield of several different mixer configurations are 
evaluated and compared with their corresponding NOx predictions. The objectives of this study 
are the following: 

1) To compare the NOx production performance for both annular and circular mixing sections 
operating at similar conditions and with similar orifice geometries. 

2) to vary the number of orifices per row in the circular mixer to minimize the NOx produc- 
tion for both round and slanted slot orifices at a given J condition, and 

3) to understand the locations for high NOx production in both annular and circular mixing 
configurations. 

Sixty different mixer configurations were analyzed with the 3-D numerical CFD code and 
the localized NOx production and overall mixer NOx production flowrate were calculated for 
each case. The NOx predictions were made by sub-dividing the overall mixing domain into 
smaller sectors. The local NOx value therefore represents the NOx production flowrate per unit 
volume (g of NOx per sec per m 3 ) since all the sectors which make the mixer domain are not 
sized equally. The mixer configurations include round holes and elongated slots with aspect 
ratio of L/W=4. Slots were modeled at both 0 degrees (aligned with the flow) and 67.5 degrees. 
The jet to mainstream momentum-flux ratio was maintained at 52 for all configurations. 


NAS A/CR— 2000-2 1 0672 


3 



Boundary conditions simulate the high speed civil transport gas turbine altitude cruise condition. 
An analysis was carried out to rank the mixer configurations with respect to both degree of 
mixing and NOx production. The mixing non-uniformity was redefined with respect to the 
stoichiometric equivalence ratio instead of the equilibrium value to investigate whether the 
mixing flowfield correlates with the NOx production. 


Mathematical Model 


3-D Flow Model 


A production 3-D combustor code, COM-3D (Bruce, Mongia, and Reynolds, 1979) is used 
that solves the turbulent reacting flow transport equations using the SIMPLE algorithm of 
Patankar and Spalding (Patankar, 1980). This program simulates turbulence by the two-equation 
k-e model (Launder and Spalding, 1974), and combustion following vaporization is determined 
by a four-step chemical reaction model based on Arrhenius and modified eddy breakup concepts. 
The transport equations for all dependent variables are of the following form as shown in Eq-1: 

dn|p r n^ - (j^jgradfc) 

where p r is the mixture density, u is the velocity, |i eff is the effective turbulent viscosity, P r is 
the effective Prandtl/Schmidt number, and S^ is the source term for the variable The 
following variables are computed by COM-3D: 1) axial, radial, and swirl velocity components; 
2) specific enthalpy and temperature; 3) turbulence kinetic energy and dissipation rate; 4) 
unbumed fuel, CO, H 2 , intermediate fuel, and composite fuel mass fractions; and 5) fuel spray 
trajectory and evaporation rate. 

The computational effort is significantly reduced by modeling a sector of the mixing section 
comprising a single orifice (or dual, generally opposed, orifices as applied to the annular mixing 
configurations). Therefore, the shape of the sector was dependent on the number of orifices 
equally spaced in the circumferential direction. Periodic boundary conditions were applied in the 
circumferential direction. No-slip and adiabatic boundary conditions were applied at the outer 
and inner walls defining the mixing section. Zero-gradient boundary conditions were applied at 
the center axis as applied to the circular mixers only. Axial gradients at the exit boundary 
condition were assumed zero. 

NOx model 



The NOx model described herein was developed by Rizk and Mongia (1993). Because NOx 
formation in the combustion liner is significantly affected by the details of the front end of the 
RQL and the subsequent admittance of air into the various downstream zones, the combustion 
liner needs to be divided into a number of regions for modeling purposes. The hybrid modeling 
technique therefore consists of using the 3-D simulation results obtained with COM-3D such as 
gas flowrate, flow averaged temperature, fuel/air ratio, in addition to the turbulence characteris- 
tics to accurately describe the flow nuances affecting the NOx production. 
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The NOx model is subdivided into three main zones: rich zone, mixing zone and lean zone. 
Each of these zones exhibit significantly different characteristics as to their role in the NOx 
production in an RQL combustion liner. Lefebvre (1984) shows that the exhaust concentration 
of NOx can be calculated in terms of mean time in the combustion zone, chemical reaction rates, 
and mixing rates. Mellor (1976) describes a model that combines a number of perfectly stirred 
reactors and uses finite-rate mechanisms to calculate NOx. Furthermore, Fletcher and Heywood 
(1971) modeled the combustor primary zone as a partially stirred reactor burning gaseous fuel 
over a distribution about the primary zone equivalence ratio, a lateral mixing reactor, and a plug 
flow reactor. Based on these studies it is clear that NOx production characteristics in an RQL 
combustor required a detailed model of the mixing flow phenomena along with detailed 
description of the chemical kinetics in the flame reaction zone. The model used in this study 
subdivides the RQL combustor into the following zones: 1) rich zone made up of two parallel 
reactors: a stoichiometric reactor to account for the combustion at near- stoichiometric conditions 
and a rich zone feed. 2) Mixer air is admitted consecutively into the quench zone to form plug 
flow reactor, and 3) dilution air is admitted into the dilution zone to form a secondary plug flow 
reactor similar to that of the quench zone. To provide an accurate estimate of the trends of 
formation and reduction of NOx, the detailed reaction model developed by Westbrook and Pitz 
(1984) was combined with extended Zeldovich mechanisms for NOx. The detailed description 
of the mixing zone model will follow. Both rich and lean zone design is categonc at this stage 
and assumptions about the modeling of these two are made as follows: 


1) The input flow field entering the quench zone is assumed to be fully reacted and in 
equilibrium. The species mole fractions for (CO, C0 2 , H ? 0, and H 2 ) were determined based 
on the given rich-zone equivalence ratio at chemical equilibrium at the prescribed operating 
conditions using JP-5 for the typical fuel properties. This assumption may not be totally 
realistic in an actual operating RQL rich section since more CO and unbumed hydrocarbon 
emissions are expected as compared to this simplified rich zone model. The additional 
emissions are also expected to affect the mixing zone reaction prompted by the mixture of the 
rich zone products with additional air. However, since the rich-zone performance is not 
within the scope of this study, the assumed simplified rich-zone model (achieving chemical 
equilibrium) is a reasonable initial condition towards analyzing the reaction in the mixing 
zone. In addition, it provides a baseline initial condition which can be applied to all 
configurations equally since it is assumed that the rich zone and mixing zone performance 
are generally coupled and the variation of one may significantly affect the performance of the 
other. The four species used in the inlet mainstream condition correspond to the four-step 
chemical reaction used in COM- 3D. 


2) No expansion is assumed in the lean zone and no area contraction is applied to the mixing 
section. The RQL liner is usually designed with an area expansion downstream of the 
quench zone in order to increase the residence time to further oxidize the CO produced in the 
rich zone. In addition, no air is added in the lean zone downstream of the quench zone. 


Based on these assumptions, the hybrid model was applied to the quench zone only. 
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Mixing zone model. 

The kinetic scheme results for the quench zone reactor show that a substantial increase in 
NOx formation occurs as the residence time in the reactor is increased. The results obtained at a 
number of pressure levels and rich zone equivalence ratios have been correlated as follows. 

f - 38,000 1 

NO q (g/kg)=L51Elx n q q e'» T < V 2 (2) 

T q is the reaction temperature in the quench reactor, in K, and Tq the residence time in the 
reactor in ms. The exponent nq is given in terms of the rich zone equivalence ratio <|> r as follows: 


nq = 0.834(J>7 0 - 36 (3) 

The hybrid model consists in grouping the results obtained from a typical COM-3D grid 
with about 50,000 to 150,000 nodes into about 2,000 larger sub-volumes comprising of several 
computational nodes representing the entire mixing section of an RQL combustion liner. By this 
means, the reacting flow model results can be accurately represented and modeled using 
significantly less computational effort. For the NOx predictions only, the model of the quench 
zone extends to x/R= 2x/H = 5 since the domain beyond the quench zone was not considered to 
be a realistic lean zone model for the given correlations. This was implemented to prevent any 
unrealistic discontinuities in the predictions near the real interface between the quench and lean 
zones (x/R= 2x/H =1) where x is the downstream axial distance from the leading edge of the 
orifice and R is the mixer radius (for circular) and H is the annulus height (for annular). 

Geometric Configuration 

Two potential design considerations for the hot section component of the high speed civil 
transport engine are the circular and annular RQL combustion liners. In a typical RQL 
combustion liner, the rich zone, quench zone, and lean zones are axially inter-connected with an 
area contraction in the mixing section to speed up mixing by increasing the mainstream 
momentum flux. The configuration details of the mixing configurations follow. The details of 
each of the 60 configurations are shown in Tables- 1 and 2. 

Circular mixer description 

The mixing section was modeled as a constant cross-sectional area throughout the axial 
direction with a single row of equally spaced orifices in the outer wall. The outer wall diameter 
is 3 inches (0.076 m) and the axial length of the mixing section extended from x/R=-1.4 to 
x/R=6. Sufficient axial distance was provided both upstream and downstream of the orifice to 
minimize the impact of the inlet and exit boundary conditions on the calculated flow structure in 
the primary domain of interest which is (0<x/R<l). The downstream limit (x/R=l) of the quench 
zone was arbitrarily defined to be x/R=l. 
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Annular mixer description 


Two different annular mixers were modeled to investigate the effect of channel height on the 
mixing flowfield and NOx production: 

a) A large channel mixer with an annulus height (H) equal to the circular mixer diameter 
where the inner and outer diameter dimensions are 13 and 19 inch respectively (0.330 m and 
0.483 m). The inner and outer walls have the same size, number, and shape orifices. The 
overall mixing domain of the computation extends from 2x/H=-3 to 2x/H=6. This 
configuration was modeled to compare the effect of jet penetration of the circular and annular 
cases with a physically the same characteristic channel height dimension. 

b) The analysis was also made with a small channel height (H) equal the circular mixer radius 
R = 1.5 inches (0.038 m). The inner and outer diameter dimensions of the mixer are 13 and 
16 inches respectively (0.330 m and 0.406 m). The inner and outer walls have the same size, 
number, and shape orifices. The overall mixing domain of the computation extends from 
2x/H=-4 to 2x/H=12. The short H annular mixer was modeled to investigate the effect of 
severe over-penetration on the mixing flowfield and NOx production. 

In order to make a valid comparison between the circular and its corresponding annular 
mixer, it became essential to first duplicate the cross-sectional area of the circular mixer modeled 
sector in the annular configuration. For this purpose, the cross-sectional area of two equal (1- 
hole) sectors in the circular geometry (corresponding to the opposing orifices) are duplicated in 
the annular configuration. The total number of orifices in the annular mixer cases were 
determined by this constraint and the physical outer and inner wall dimensions mentioned earlier. 

Note that based on the constraints imposed above, both the circular and annular cases with 
the same orifice shapes will yield the same hole size for a constant operating condition. Orifice 
Cd is assumed to be 0.7 throughout the entire study. The description of the orifice types and 
mixer configurations for all 60 configurations analyzed herein are as shown in Tables- 1 and 2. 
Note that the actual number of holes per row are detailed separately from the equivalent number 
of holes. The equivalent number of holes is used to compare both annular and circular mixing 
configurations and is defined as the number of orifices per row per equivalent sector cross- 
sectional area as shown in Eq-4. 


W„ = - 


N c ( Do?-DC) 
2(Do e 2 ) 


(4) 


where, N^ = the equivalent number of annular orifices, N c is the number of orifices in a circular 
can, Do a and Di a are the outer and inner annulus radii of the annular mixing configuration, and 
Do c is the outer diameter of the circular mixer. 


The equivalent sector is that given by the cross-sectional area of 2 opposing orifices in the 
circular mixer. In the circular mixer, the actual and equivalent number of holes is the same, 
however, for the annular mixer, the equivalent number of holes is significantly reduced since the 
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overall area of the annulus is larger than that of the circular mixer. In other words, the equivalent 
number of holes indicates how many orifices in the annular configurations will it take to equal 2 
opposing orifices in the circular mixer. Also note that the calculated equivalent number of 
orifices is per wall, recognizing that both inner and outer walls in the annular configuration have 
the same number of orifices. The equivalent number of orifices for all configurations are shown 
in Tables- 1 and 2. Note that in all of the annular configurations except for three different 
configurations, the orifice arrangement for the annular mixer are directly opposing orifices. The 
orifices in both the inner and outer walls for HO-34, 40, and 41 are in a staggered arrangement. 

The computational domain was discretized into 50,000 to 150,000 finite control volumes 
arranged with about 70 to 120 nodes in the axial direction, and 30 to 40 nodes in both the radial 
and tangential directions. The grid was denser near the orifice to resolve the high velocity and 
temperature and equivalence ratio gradients resulting from the inlet of the crossflow jet. An 
orthogonal view of a typical grid arrangements for both circular and annular mixer types are 
shown in Figures- 1 and 2. The grid is configured to allow smooth progressive volume change 
between adjacent control volumes to help speed up the convergence of the solution. In the 3-D 
numerical model, a secondary grid (staggered grid) is interpolated from that shown in Figures- 1 
and 2 to obtain the boundaries upon which the vector quantities are acting. In other words, the 
staggered grid becomes the control surfaces of the micro control volumes that constitute the 
inner volume of the mixing section. 

A total of 60 circular and slanted slot orifice configurations were analyzed as shown in 
Tables- 1 and 2. For the circular configurations, the blockage is defined as the circumferential 
projection of the orifice divided by the spacing between the orifice centers. The T.E. x/R is the 
non-dimensional axial location of the trailing edge of the orifice with respect to its leading edge 
(x/R=0). For the annular configurations, the blockage is defined as the circumferential 
projection of the orifice divided by the spacing between the orifice centers evaluated at the mid- 
area radius. The mid-area radius is the radius location that equally splits the cross-sectional area 
in 2. The orifice trailing edge non-dimensional parameter as applied to the annular configuration 
is T.E. 2x/H where the H is the annulus height. Both the blockage and the trailing edge 
nondimensional parameters are shown in Tables- 1 and 2. 

In this study, the following parameters were kept constant throughout this analysis. The 
control of these parameters is essential for a valid comparison of both mixing non-uniformity and 
NOx production between the appropriate configurations. 

1) rich zone equivalence ratio = 1.80 

2) lean zone equivalence ratio = 0.416 

3) overall mixer pressure = 203psia (14.1 atm) 

4) jet temperature = 17 10 R (950 K) 

5) mainstream temperature = 3895 R (2164 K), from the chemical equilibrium code (CEC) 
developed by NASA 

6) mixer diameter = 3 inch (0.076 m) (circular) 

7) Di=13 inch (0.330 m) and Do = 16 inch (0.406 m) for the small annular mixer 

8) Di=13 inch (0.330 m) and Do = 19 inch (0.482 m) for the large annular mixer 

9) jet to mainstream mass flow rate MR = 2.96 

10) jet to mainstream density ratio DR = 2.28 
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11) jet to mainstream momentum-flux ratio (J) = 52 

12) overall mixer total pressure loss = 2.69% 

The following were allowed to vary: 

1) orifice blockage and T.E. x/R values depending on the orifice shape and orientation 

2) number of orifices per row 

3) orifice shape 


Modeling specifications 

The species mole fractions, temperature, and velocity profiles exiting the rich zone and 
entering the mixer were assumed to be uniform across the inlet cross-section of the mixing 
section. The air jet flow was characterized by a radial, uniform flow across the orifice effective 
area. COM-3D is not a body conforming code therefore accurate modeling of the orifice shape 
was obtained by a stair-stepping approximation using 80 to 300 control surfaces. This number of 
nodes/orifice allowed accurate modeling of both round and slot geometries. The assumption of 
uniform mass injection/area for the orifice air entry is applied in the mathematical model in all 
the configurations. The turbulence kinetic energy of the mainstream and jet flows were 0.3% of 
the square of the mean velocities. The turbulence length scales of the mainstream flow were 2% 
of the mixer diameter, and the turbulence length scale of the jet was of the order of the orifice 
diameter. The results from COM-3D were post processed into about 2000 sub-volumes, about 
20 in the axial direction, 10 in the radial, and 10 in the circumferential direction. This input was 
used to generate the NOx predictions based on the flow characteristics. 

A typical numerical solution took about 500 iterations for full reacting flow convergence 
with overall mass flow residuals of 0.05% of the total mixing section mass flowrate. All 
solutions were obtained using the Cray C-90 and a typically converged solution took about 2 hrs 
of CPU time. 

Result representation 

The temperature, local NOx, and normalized equivalence ratio plots are also shown in the 
figures in Appendix-A for all 60 configurations described in this report The 3-D perspective of 
these plots places the observer downstream of the mixer looking upstream at the sector of the 
mixer. Note that configurations with many orifices per row tend to produce very thin single- 
orifice sectors. When this occurs, several sectors are plotted adjacent to each other in order to 
better illustrate the flow field (see Figure-A58). Five different radial-tangential planes have been 
chosen at x/R = 2x/H = 0.08, 0.25, 0.5, 0.75, and 1.0 showing the flow development and NOx 
production, mixing flowfield, and reaction temperature throughout the mixing section. The axial 
direction is plotted vertically and the flow moves from bottom towards the top. 
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To the left of each plot in Appendix-A is the temperature development in K. Nine different 
contours are represented in this plot ranging from 900K to 2500K. 

The middle plot in Appendix-A figures is the localized NOx production in log(gr of NOx / 
{sec m 3 }). The contour scale ranges from -4 to 4 and the shades indicate that the red contours 
are locations of the highest NOx production and consequently the blue indicates the lowest NOx 
production regions. 

The parameter (f) defined in Eq-5, applied to the equivalence ratio scalar quantity ($), is 
plotted to the far right in all figures in Appendix-A. Note that the equivalence ratio is a 
conserved scalar. The value of (f) varies from 0 to 1, where 0 (blue) is the value of the unmixed 
jet equivalence ratio and 1 (red) is the value of the mainstream flow equivalence ratio. Note that f 
= 1-0, where 0 is as defined previously (Holdeman, 1993) and used elsewhere also. Also note 
that the definition of (f) does not apply to the temperature distribution in the reacting flow 
solutions since the temperature in some areas of the mixing region rise above that of the 
incoming mainstream flow due to the chemical reaction of the rich-zone exit composition. The 
equilibrium (f) value was calculated in a purely adiabatic system at any location downstream of 
the jet injection. 


/= 


§ main ^ jet 


(5) 


Results and Discussion 

The mixing performance for all configurations analyzed in this study were ultimately 
evaluated at x/R = 2x/H =1. The flow mixing non-uniformity throughout the mixer volume (0 < 
x/R = 2x/H < 1) is also quantified since complex shear structures are present in the flow field 
near the jet entry region. 

Two different methods were used to analyze the numerical results: observational analysis 
and statistical analysis. Both of these are described in detail below. 

Observational analysis 

Temperature, equivalence ratio, and NOx production field solutions were analyzed 
qualitatively using the figures shown in Appendix-A. Mixing non-uniformity is clearly shown 
by the variation in the color distribution in both temperature and equivalence ratio distribution 
contours. The temperature plots indicate locations, especially behind the orifice, where the 
reaction between the mainstream and jet flow produce the expected temperature rise. The NOx 
plots in Appendix-A show the critical NOx production locations. Note that the normalized 
equilibrium equivalence ratio for all configurations is 0.253. Similarly, the equilibrium 
temperature for all configurations is 1598K. 
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Statistical analysis 


The statistical analysis involves the detailed characterization of the mixing flowfield using 
several different techniques. These techniques characterize the mixing non-uniformity only. 
Three different statistical methods are reported in this report and are shown in the figures in 
Appendix-B: 

a) The performance of the mixing section at the radial-tangential planes at x/R=l and 2x/H 
were quantified using an area weighted planar deviation parameter (AMIX) given by Eq- 
6. The deviation used in AMIX is calculated relative to the stoichiometric value of f 
which equals 0.556 for all configurations. The stoichiometric value is used as a 
reference because most NOx is produced in the rich regions of the flame. A large AMIX 
value indicates a more uniformly distributed flow and lower NOx production tendency. 
It is also important to note that AMIX does not correct for the bias introduced in the 
region of the mixer where the air is being injected through the orifice. However, AMIX 
is only evaluated at the exit of the mixer x/R= 2x/H =1 and this plane is downstream of 
the air injection in almost all configurations summarized in this report Exceptions to 
this are configurations HO-19, 59, and 60 which have a T.E. 2x/H larger than 1.0 (see 
Tables-1 and 2). The area-weighted non-uniformity results are shown in Tables-3 and 4 
in configuration nomenclature order and in Tables-5 and 6 sorted in order of mixer type. 
The AMIX is shown as a function of mixer axial location in the top plot in the figures in 
Appendix-B. 


AMIX = 


1 V 1 , f §jk §stoich 1 1 

vf J 


( 6 ) 


b) The mass flow weighted planar deviation parameter MMIX (also evaluated at x/R— 1 
or 2x/H) as defined in Eq-7 is also used to evaluate the mixing region. This parameter is 
similar to the area weighted parameter in Eq-6, with the added density and velocity 
weighting terms. MMIX parameters for all configurations are shown in Tables-3, 4, 5 
and 6. As with the AMIX definition, large MMIX values correspond to better mixing 
and lower NOx production tendency. Axial profiles of MMIX are shown in the top plot 
in the figures in Appendix-B. 


2 


MMIX = 
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(7) 


c) The flow field was also evaluated by calculating a numerical volume integration NVI 
parameter for each mixing section (0 < x/R < 1 or 0 < 2x/H < 1) as shown in Eq-8. 
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^Volume ijk 


NVI = Volume Fraction! bin = 


A_ 
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J 'Volume# 

fjet 


( 8 ) 


The NVI parameter allows analysis of the entire flow field and is more descriptive of the 
overall flow phenomena in the mixing section as compared to a planar deviation analysis 
shown in Eqs-6 and 7. This NVI analysis was only performed on the equivalence ratio 
distribution. The entire range of the normalized parameter f (from 0 to 1) was sub- 
divided into 200 equal size bins and the volume of the computational control volumes 
corresponding to the value of (f) at a certain bin size (fj to fi+^f) was integrated as 
shown in Figure-3. The integrated volume in each bin was normalized based on the 
entire analyzed mixer volume, thus obtaining the normalized volume fraction. The 
volume fraction for each bin was plotted as a histogram as a function of f. In addition, 
the incremental partial volume per bin was integrated from f = 0 to 1 and the cumulative 
volume fraction was obtained for all the 60 analyzed configurations as shown in Figures 
4 and 5. If the histogram is integrated from f=0 to f=l, the value of the cumulative 
volume fraction will be 1 since the overall volume itself is also normalized to 1. 

The shape of the volume fraction histogram was also characterized by the definition of 
B(+) and B(-). The value of B(+) is the "area determined" distribution half width above 
fstcich such that the integrated area under the histogram above f stoich is 1/2 that of the 
overall histogram area above f stoich . The same applies to B(-) for the area below f stoich . 
The definitions of B(+) and B(-) are shown pictorially in Figure-3. The lowest NOx 
tendency yields highest B(+) and B(-) indicating that most of the flowfield within the 
mixing volume is either too rich or too lean to produce significant amount of NOx. The 
histogram shape value SUM-B is defined as the sum of B(+) and B(-). The shape of the 
volume histograms for all 60 configurations are shown in the bottom plot in the figures 
in Appendix-B. 

Effect of the change in the number of orifices/row on jet penetration 

In this comparison it is important to reiterate that the comparison is made without changing 
the overall ACd of the row of orifices and that the variation of the number was accomplished by 
changing the size of the orifice to keep the same row ACd. This was necessary in order to keep 
the jet-to-mainstream mass-flow ratio constant for all configurations. In addition, the shape of 
the orifices, whether round or slot, was preserved and the hole became smaller for a larger 
number or orifices. It is also important to note that for certain annular configurations, a 
maximum number of orifices is associated with the case and this means that the web between the 
orifices in the circumferential direction becomes zero. 

Jet penetration as shown in Tables-3, 4, 5 and 6 is defined from the outer wall of the mixer 
in both annular and circular mixers. The jet penetration is expressed as a non-dimensional 
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quantity r/R for the circular configurations and 2r/H for the annular configurations. The jet 
penetration is also calculated by plotting the normalized equivalence ratio as a function of radius 
along the axial-radial plane through the center of the orifice. The jet penetration is therefore 
evaluated at x/R=l (circular) or 2x/H (annular) at the radial location where the value of f is 
minim um thereby indicating the location of the center core of the jet. 

The variation in the jet penetration for both annular and circular mixers are shown in 
Figures-4 and 5 respectively. The abscissa value in both plots is the equivalent number of 
orifices which has been defined previously in Eq-4. These figures also show the difference in 
the jet penetration as a function of orifice shape as a parameter. In addition, Figure-5 indicates 
the effect of the change in annulus height on the jet penetration. . The results indicate the jets 
have penetrated to the center of the mixer in almost all cases of the small annulus height mixers 
with the round orifices (labeled as 16ANN-RO, which indicates 16 inch Do, annular mixer with 
round holes). The 16ANN-SS results however indicate moderate decrease in jet penetration with 
the increase in number of orifices. Note that the SS indicates the 67.5° slanted slot L/W=4 
orifice type. The large annulus mixer 19 ANN- AS and 19ANN-RO results for the aligned slot 
L/W=4 and round orifices respectively indicate also moderate decrease in jet penetration with 
increasing number of orifices. Note that the AS nomenclature is used for the flow aligned slot 
L/W=4. The trends shown in Figure-4 are expected, especially for the difference between the 
slanted slot and round holes for the same mixer type. The slanted slot has been shown 
previously by Oechsle, Mongia, and Holdeman, 1992, 1993, and 1994, to have lower jet 
penetration due to the transfer of the radial momentum to circumferential recirculation as the jet 
penetrates into a circular mixer. The same phenomena is shown herein for annular mixers as 
well. Note that the jet penetration for the staggered hole arrangement are not plotted in Figure-4 
since the definition of the jet penetration does not hold when the jets penetrate beyond the mid- 
plane of the annular mixer. 

The results in Figure-5 for the circular mixers also show the expected trends for the AS, RO, 
and SS orifices (flow aligned slot L/W=4, round, and 67.5° slanted slot L/W=4) respectively. 
The decrease in jet penetration with increase in orifice number is weaker for the strongest 
penetrating configuration (the aligned slot) as compared to the weakest (the slanted slot). These 
results are confirmed by other studies in the literature. It is however important to note that the 
orifices in the circular mixer have much more difficulty penetrating to the center core as 
compared to the annular cases. For the cases where the jets in the circular mixer penetrate to the 
center core (which is a singular axis in the flow field), the mainstream flow produces great 
resistance. This however is not shown in the annular mixers since there is no singularity at the 
center of the mixer and therefore the jets can more readily penetrate completely to the center as 
shown in Figure-4 with values of jet penetration equals 1. 

Effect of the change in number of orifices/row on AMIX and MMIX 

The AMIX variations as a function of the equivalent number of orifices per row are shown 
in Figures-6 and 7 for the annular and circular mixers respectively. Reiterating that the large 
value of AMIX has a lower NOx production potential, the results in Figure-6 show that the 
orifices with high jet penetration such as the 16ANN-RO, and the 16ANN-RO-STAG have 
minimum AMIX at 22 eq holes/row. However, orifices with lower jet penetration such as the 
19ANN-RO, 19ANN-SS, and 19ANN-AS minimized AMIX at much lower eq number of holes 
(about 8). This same trend appears to be similar in Figure-7 results for the annular mixer type. 
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The orifices producing stronger penetration seem to have minimized AMIX at higher number of 
orifices per row (8, 12, and 18) for the SS, RO, AS respectively. 

The MMIX variations as a function of the equivalent number of orifices per row are shown 
in Figures-8 and 9 for the annular and circular mixers respectively. The results appear to be 
similar to the AMIX results. Note that the a small MMIX indicate a greater tendency to produce 
NOx in the mixing section. The orifices that penetrate to the center core as the number of holes 
is decreased tend to have a lower MMIX at a higher number of equivalent holes per row (for the 
cases shown herein, about 20 to 30 eq holes/row). Configurations with shallow penetration and 
fewer holes have higher values of MMIX (see Figures-8 and 9). It is also evident from the 
Figure-10 that the jet penetration for minimum AMIX value is about 0.7 to 0.8 for the annular 
mixers while being about 0.6 or less for the circular mixers. The correlating trends between the 
MMIX parameter and the jet penetration (shown in Figure- 11) are also similar to that of the 
AMIX case for both the annular and circular mixers. 

Effect of the change in number of orifices/row on SUM-B 

SUM-B variations as a function of the equivalent number of orifices per row are shown in 
Figures- 12 and 13 for the annular and circular mixers respectively for mixer domain of 0<x/R<l 
or 0<2x/H<l. It is important to reiterate that the largest values of SUM-B tend to produce less 
NOx. The results indicate that in both annular and circular mixers, as the number of orifices is 
increased, the value of SUM-B also increases for all orifice shapes. These results are expected 
since, although the values are normalized for the unit volume in the mixing section, as the 
volume sample decreases, the likelihood to encounter larger gradients in the flow also increases. 
This also indicates that mixing uniformity degrades as the volume decreases which appears to be 
common for all configurations with large number of orifices per row. These results support the 
results shown in Oechsle and Holdeman (1995) that volume is essential for adequate mixing. 
Therefore it is expected that the SUM-B will increase for a reduced number of orifices. The 
results shown in Figure- 14 indicate good correlation between the SUM-B parameter and the jet 
penetration for the annular and circular mixers. Decreasing SUM-B appears to be prevalent for 
all configurations exhibiting stronger penetration. Note that the 16ANN-RO data values in 
Figure 14a appear not to follow this trend. This is expected since all these cases are over- 
penetrating to the core of the mixer. The results in Figure-14 indicate that as the jet penetrates to 
the core, a larger fraction of the flow field with a fuel concentration approaching the stoichiomet- 
ric value increases, thereby resulting in a tendency to produce more NOx. It is important to note 
that although the volume fraction histogram defines the partial volume in the mixer with a 
specific f value near the stoichiometric value, this parameter is not a function of the residence 
time which is an important factor in the NOx production. Therefore the correlation between 
SUM-B and NOx reduction should carry the assumption that the flow field mean velocity is 
generally constant in the axial direction. This, of course, may not be applicable in some cases, 
especially in flow fields with large recirculation vortices where pockets of the mixture appear to 
reside longer at a certain value of f due to the transient time within the recirculation eddy. 

Effect of the change in number of orifices/row on mixer NOx production 

In this comparison, the NOx cumulative concentration is calculated up to the exit boundary 
of the mixing section. As mentioned previously, the exit boundary for the circular mixer is 
x/R=l and for the annular mixer is 2x/H. The NOx concentration (CUM NOx) in ppm is 
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accumulated throughout the mixing section, including any NOx formed upstream of the x/R - 
2x/H=0 inlet boundary. This is important since the jet flow in many of the over-penetrating 
mixing configurations tend to recirculate upstream of the x=0 boundary mixer and mix with the 
mainstream flow. When this occurs, temperature rise is expected in the region upstream of die 
designated mixer inlet and thus NOx is formed. It is also important to note that although a 
designated exit boundary has been imposed in this study, additional mixing and reaction can 
occur beyond this limit, therefore producing NOx that is not accounted in the CUM-NOx value 
described above. Therefore, a second NOx parameter is defined (Global EINOx) which is the 
NOx emission index for the entire computational domain. This parameter is considered valuable 
since it accounts for the NOx that is produced downstream of the mixer up to a value of x/R = 
2x/H = 6. Note that the sample length of the mixer is the same for both the circular and annular 
mixers. This is essential since NOx is significantly dependent on the residence time. It is also 
assumed that this length is sufficient to allow full reaction of the available dissociated Oxygen 
atoms enough to assume that no NOx will be produced further downstream. It is therefore 
believed that the global EINOx is a more descriptive value of the NOx production in a mixing 
region. 

The effect of the variation of the number of holes on the CUMNOx is plotted in Figures- 15 
and 16 for the annular and circular mixers respectively. The results indicate no simple 
relationship between the NOx production in the mixer as the number of equivalent holes in 
increased in the range analyzed. It appears that the NOx production is somehow dependent on 
the jet penetration. The results shown previously indicated increase in the number of orifices 
resulted in a decrease in jet penetration for constant MR, J, and DR. It is also recognized that 
once the jet penetration has reached the core, further reduction in the number of orifices could 
have secondary effects on the mixing and NOx production. The results shown in Figure- 17 
indicate that there are multiple values of the jet penetration where NOx production is low. It is 
therefore apparent that there is no unique correlation between the NOx production in the mixing 
section (CUM NOx) with jet penetration. However, differences in the behavior between circular 
and annular mixers as shown in Figure- 17a and b indicate that possible change in the trends from 
over-penetration to under-penetration could possibly cause different NOx production trends as 
well. 

The effect of the variation of the number of holes on the Global EINOx value is plotted in 
Figures- 18 and 19 for the annular and circular mixers respectively. Although there does not 
appear to be a simple relation between the overall NOx production and increasing number of 
equivalent orifices, there is however a possible qualitative explanation for the trend. All of the 
orifices in all of the mixer types indicate the same general trend in the NOx production as the 
number of holes is changed. These results indicate that jet penetration is definitely a strong 
factor in the NOx production. If the jet penetrates too much, the jet flow will be assimilated into 
the mixer flow in the center leaving a largely poorly unmixed region behind the orifice where 
NOx is produced. Similarly, if the jet penetration is low, a poorly mixed region will be created 
near the core of the mixer. To understand the trends shown in Figures 18 and 19, the observer 
must follow the trend from the under-penetrated side (moving in the abscissa from the right to 
left) as the number of equivalent orifices per row is continuously decreased. The trends will be 
described in detail below in three phases: 

Phase-a) As the number of orifices is decreased, the jet penetration increases (see Figures-4 

and 5) slowly towards the optimum level and, the NOx production decreases to its lowest 
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value. The cases that show this trend are the 19 ANN- AS at 18 equivalent holes per row in 
Figure- 18 and all the trends shown in Figure- 19 for the circular mixer at 30, 20 and 10 
holes/row for the AS, RO, and SS hole shapes respectively. It is interesting to note that the 
minimum NOx is attained for these configurations at jet penetration between 0.4 and 0.55. 

Phase-b) When optimum jet penetration is reached, a further decrease in the number of 
orifices per row increases the jet penetration and resulting in increased NOx. Cases GANN- 
AS clearly indicate this trend in Figure-18 and for the CAN-AS and CAN-RO (ranges from 30 
to 16 holes/row) and (20 to 10 holes/row) respectively shown in Figure-19. The NOx 
production in phase-b generally increases as the jet penetration increases beyond the optimum 
value. Note that the NOx production peaks in this phase, however, this peak may or may not 
correlate with the jet over-penetration to the center core as shown in Figures-4 and 5. 
However, it is evident that further increase in the jet penetration beyond the high NOx peak in 
phase-b, actually lowers the NOx production. 

Phase-c) This is the final phase of the trend. When a reduction in the number of orifices 
produces further increase in jet penetration, two things may happen: 1) the jet may have 
already penetrated to the center core in which case it cannot go any further, or 2) the jet is still 
approaching the center core. It is apparent that with either of these cases, the increase in the 
jet penetration produces a decrease in NOx formation. The decrease in NOx formation is 
primarily aided by secondary mixing phenomena that has been previously neglected due to the 
small effect it has during phases a and b. For example, the annular configuration benefits 
from this over-penetration by forming vortices in the circumferential direction caused by the 
collision of the opposing jets. However, in the circular mixer, this can not occur since the 
volume in the center core of the mixer is a singular point. Nevertheless, further mixing is 
aided in the circular mixer, especially in the region behind the orifice. In this region, the 
upstream recirculating flow created by the colliding jets re-enters the mixer usually behind the 
orifice since the jet flow is occupying the center core. The effect of the grossly over- 
penetrated cases is shown schematically in Figure- 20. 

The illustration of the three phases in two different orifice shapes and mixer types are shown 
in Figure-2 la. The other configurations not mentioned in these three phases only exhibit part of 
the full trend just described. The reason for this is that in many cases, the maximum number of 
orifices was reached and no further orifices could possibly be fit in the annular configurations. 
This was the case with 16ANN-SS, 16ANN-RO, 19ANN-RO. The other cases were not 
investigated due to time constraints. Also it is important to note that the value of the number of 
orifices for minimized NOx formation is different for different orifice shapes and mixer types. 
This means that the number of required orifices to produced the lowest NOx is generally higher 
for a round hole or aligned slot which generally have stronger jet penetration. In similar manner, 
weaker penetration orifice configuration such as the slanted slot will require fewer number of 
holes to obtain an optimum design. 

It is difficult to access whether optimum number of orifices was found for all the orifice 
configurations, especially since not all three phases of the trend are shown for some of the cases. 
However configurations such as the 19ANN-AS, 16ANN-SS, and all the circular mixer orifices 
indicate that generally the lowest NOx attainable form either the annular or circular mixers is 
about 2 NOx El and the overall variation of NOx is about 7 times the minimum. Therefore an 
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optimization study is recommended for any RQL mixer design since it has the potential of 
having a significant NOx reduction at a given operating condition. 

The plots shown in Figure-2 la also indicate that the comparison between the annular and 
circular mixers indicate similar trends as mentioned above. All three phases of the trend are 
shown for the aligned slot L/W=4. However, the trends are shifted with respect to each other. 
The annular configuration shows the trend shifted towards lower equivalent number of orifices 
per row as compared to the circular results. It is important to note that in this comparison the 
physical dimension between the opposing orifices is the same (3 inches) for both the circular and 
annular mixers. It is also anticipated that in this comparison (where J, MR, and DR are 
constant), similar jet penetrations are expected for the same orifice shapes. However, the jet 
penetration into the circular geometry is inherently different compared to that of the annular case. 
This is probably the main reason for the shift in the trends in Figure-2 la. Also note that the 
lowest NOx was over 50% less in the annular case as opposed to the circular mixer for a given 
orifice shape. 

If the comparison is made for a given mixer geometry type with different orifice shapes as 
shown in Figure-2 lb, similar trends are found. The difference in the jet penetration for the round 
orifices as opposed to the slanted slot are significant (see Figure-4) for the annular configuration. 
However, the trends shown in this figure are not complete as the ones shown in Figure-2 la. 
Most of the trend for the 16ANN-RO (16 inch Do annular configuration with the round orifice) 
over-penetrated throughout the entire analyzed range from 4 to 30 equivalent number of orifices, 
therefore the portion of the trend shown is the phase-c part with a small part of phase-b. The 
16ANN-SS however, due to the shallower jet penetration, exhibits the phase-a from 12 to 20 eq 
holes/row, phase-b between 8 to 12, and phase-c below 8 eq holes/row. The significant shift in 
the trends for the round to the slanted slot are again explained by the difference in the jet 
penetration as a function of J and MR. Note that for a given mixer, the variation of orifice shape 
(round to slanted slot, see Figure-2 lb) appears not to have significantly affected the minimum 
NOx value (EI=4). This means that if the designer has the freedom to vary the number of 
orifices per row in a given design, any orifice shape will probably be adequate in obtaining an 
optimum low NOx value. 

Similar trends are shown for the different orifices in the circular mixer in Figure- 19. Both 
the round hole and the flow aligned L/W=4 slot have complete trends exhibiting the three phases 
in the shown spectrum. The 67.5° slanted slot L/W=4, however, has the lowest jet penetration of 
all three orifice shapes plotted in Figure- 19 and therefore has an incomplete trend with part of 
phase-c missing. Note that the shift in the trends for all three orifices moves to the left for 
orifices with lower jet penetration as a function of J and MR. 

Correlation between the mixing parameters and NOx production 

The results between the correlation of the global NOx El and the different mixing non- 
uniformity parameters are shown in Figures- 22, 23, and 24 for the phase-b results only. The 
phase-b results were chosen for several reasons, including the fact that most of the data points 
are within this regime. It is believed that all three phases will have different correlation factors, 
therefore only one is chosen in this comparison. The results show definite trends in the NOx 
formation with change in the mixing non-uniformity as expected. Note that in all plots in 
Figures- 22 through 24, the NOx production decreases with increasing mixing non-uniformity 
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parameters (AMIX, MMIX, and SUM-B). As the AMIX, MMIX, and SUM-B parameter 
increase, the flow field in the mixing region is more stratified away from the stoichiometric value 
and therefore has less tendency to produce NOx. The best correlation was found to be the SUM- 
B value. This is not surprising since this parameter considers the NOx production throughout the 
entire mixer as opposed to the AMIX and MMIX values which only consider the region 
downstream of the jet injection. Also note that Figure-22, 23 and 24 indicate that the orifice 
shape is a variable to be considered in the correlation since it significantly affects the jet 
penetration. 


Conclusions 

1) The jet penetration is significantly dependent on the mixer geometry whether a circular or 
annular configuration. The mixing flowfield and NOx production are different in the circular 
and annular mixers. Significant NOx production was shown in the region behind the orifices in 
the circular and annular mixers at locations beyond x/R=l for configurations with deep jet 
penetration. The cases with shallow penetration tend to produce NOx near the center core of the 
mixer. Optimum jet penetration minimizes the NOx production in these two regions. NOx 
production in the shear layer between the jet and the mainstream flow will generally be present 
and difficult to eliminate. 

2) The statistical parameters AMIX, MMIX, and SUM-B appear to correlate better with the NOx 
production when the stoichiometric reference value is used instead of the equilibrium value. 

3) The NOx production variation as a function of the equivalent number of orifices was found to 
have a common trend development for all orifice configurations and mixer types. The trend has 
been divided into three phases which cover the range between over-penetration, optimal 
penetration, and under-penetration. Each phase has an associated NOx trend as a function of the 
variation of number of orifices. Minimum NOx formation was found in the interface between 
phase-a and phase-b which coincides with the optimum jet penetration. Optimum penetration 
was between 40 to 60% of the characteristic dimension between the outer wall and the center of 
the mixer. The variation of this optimum penetration was surprisingly low between different 
orifice configurations and mixer types. 

4) The comparison of the performance of the aligned slot L/W=4 in both circular and annular 
mixing configurations for constant J, MR, DR, and orifice ACd/row indicate that the annular 
mixers produce significantly lower NOx as compared to the circular counterpart. However no 
measurable difference was found in the lowest NOx production for the round orifice and the 
67.5°slanted slot L/W=4 when both are used in an annular mixer. This also indicates that if the 
designer has the freedom to vary the number of orifices/row in a given mixer design, any orifice 
shape will probably be adequate to obtain a low NOx production value. This statement is 
however qualified by the fact that the jet penetration is crucial and that there is a minimum 
number of orifices tare required to achieve the optimum jet penetration. 
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Table-1. Configuration description for cases HO-1 through HO-30 
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Table-2. Configuration description for cases HO-31 through HO-60 
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Table-3. Mixing uniformity and NOx production results for the cases HO-1 through 30 
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0339 

0.600 

52 

14 

CIR 

alg slot 

0.87 

0.181 

0364 

0345 

0333 

0.335 

65 

10.7 

3 

0.67 

0382 

0.507 

53 

18 

CIR 

alg slot 

0.803 

0.277 

0.351 

0.628 

0324 

0344 

49.7 

8.8 

3 

0.6 

0321 

0.448 

54 

22 

CIR 

alg slot 

0.794 

0.41 

0.364 

0.774 

0334 

0367 

49.1 

8.07 

3 

034 

0355 

0.405 

55 

30 

CIR 

alg slot 

0.788 

0.423 

0.392 

0.815 

0364 

0399 

38.5 

4.76 

3 

0.46 

0.414 

0347 

56 

34 

CIR 

alg slot 

0.791 

0.432 

0.403 

0.835 

0372 

0.404 

35.8 

6.33 

3 

0.44 

0.440 

0325 

57 

38 

CIR 

alg slot 

0.795 

0.413 

0.404 

0317 

0375 

0.41 

40.1 

10.1 

3 

0.44 

0.468 

0309 

58 

28 

AN 

alg slot 

0.756 

0352 

0306 

0.858 

0.453 

0.43 

40.9 

10.6 

19 

0.44 

0.792 

0360 

59 

2 

AN 

alg slot 

0.765 

0.243 

0399 

0.642 

0398 

0.444 

49.2 

531 

19 

1 

0313 

1355 

60 

1 

AN 

alg slot 

0.724 

0.216 

0368 

0384 

0379 

0.443 

34.93 

5.57 

19 

a. /-A 

1 

0.147 

1372 


Table-4. Mixing uniformity and NOx production results for the cases HO-31 through 60 
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HO ; 

fteq 1 

Mixer 

Hole < 

Cum Vol 1 

B(+) 1 

B(-) ! 

5UM-B , 

AMIX 1 

V1MIX ( 

Cumul 1 

EINOx 

Do 1 

Pen r/R 1 

Blockage ' 

T.E.X/R 

' 

holes 

Type i 

type ( 

@fstoi 





i 

NOx, ppm 

Global 

inch 1 

From Do 

—1 

T.E. 2x/H 

60 

1 

AN 

alg slot 

0.724 

0.216 

0368 

0384 

0.379 

mm\ 

34.93 

537 

19 

wmmM i 


1372 

59 

2 

AN 

alg slot 


0.243 

0399 

0.642 

0.398 

0.444 

49.2 

5.51 

19 

i 

0.213 


43 

6 

AN 

alg slot 

0.706 

0.339 

0.413 

0.752 

0347 

0394 

47.5 

9.6 

19 

i 

0366 

WESSami 

47 

10 

AN 

alg slot 

0.732 

0.432 

039 

0.822 

0347 

0373 

46.2 

53 

19 

0.76 


mmmi 

48 

14 

AN 

alg slot 

0.748 

0.379 

EX29I 

EESII 

0.425 


16.8 

2.52 

19 

■HI 

0.561 


49 

18 

AN j 

alg slot 

0.756 

0.331 

0.51 

0341 

0.453 

0.423 | 

14.6 

1.6 

19 

■911 

0.634 

0.448 

50 

22 

AN 

alg slot 

0.758 

0.326 

0.5 

0.826 

E2ZEII 


9.9 

2.4 

19 

039 | 

0.701 

0.405 

45 

26 

AN 

alg slot 

0.768 

0.368 

0305 

0373 

0.455 

0.438 

14.7 

3.4 

19 

0.4 

0.763 

0373 

58 

28 

AN 

alg slot 

0.756 

0.352 

0306 

0.858 

0.453 

0.43 

40.9 

10.6 

19 

0.44 

0.792 

0.360 

19 

4 

AN 

round 

0.799 

0.15 

0.412 

0362 

0399 

0.453 

26.8 

63 

16 

1 

0333 

1.052 

11 

6 

AN 

round 

0.789 

0.181 

0.43 

0.611 

0356 

0.404 

243 

4 

16 

1 

0.405 

0352 

12 

8 

AN 

round 

0.772 

0.247 

0.426 

0.673 

0352 

0.4 

26.9 

53 

16 

1 

0.465 

0.735 

13 

10 

AN 

round 

0.762 

0.269 

0.421 

0.69 

0345 

0393 

28.4 

5.1 

16 

1 

0324 

0.662 

14 

12 

AN 

round 

0.782 

0.229 

0.413 

0.642 

0342 

0392 

33.2 

5.6 

16 

1 

0372 

0.602 

15 

14 

AN 

round 

0.759 

0.264 

0.412 

0.676 

034 

0391 

33.5 

53 

16 

1 

0.616 ' 

0.556 

25 

16 

AN 

round 

0.767 

033 

0.404 

0.734 

0336 

0388 

35.6 1 

63 

16 

1 

0.662 

0323 

26 

18 

AN 

round 

0.77 

0.4 

0356 

0.756 

0332 

0384 

33.8 

7.1 

16 

1 

0.700 | 

0.492 

27 

20 

AN 

round 

0.77 

0.365 

0386 

0.751 

032 

0377 

15.1 

12 

16 

1 

0.737 

0.466 

28 

22 

AN 

round 

0.78 

0.357 

0386 

0.743 

032 

0376 

15.1 

9.2 

16 

1 

0.776 

0.446 

29 

24 

AN 

round 

0.784 

0.357 

0386 

0.743 

0324 

0373 

15.4 

9.2 

16 

1 

0309 

0.426 

33 

30 

AN 

round 

0.777 

0.352 

0.43 

0.782 

0361 

0368 

15.5 

7.1 

16 

0.91 

0.904 

0381 

40 

2 

AN 

round 

0.802 

0.233 

0.43 

0.663 

035 

0395 

36 

73 

19 

1 

0.468 

0.744 

39 

6 

AN 

round 

0.794' 

0.343 

0395 

0.738 

0337 

0375 

62.8 

14.4 

19 

035 

0304 

0.426 

20 

6 

AN 

sit slot 

0.663 

0.3 

0.417 

0.717 

0346 

0.39 

17.6 

7 

16 

0.95 

0.681 

0394 

21 

8 

AN 

sit slot 

0.647 

0.422 

039 

0.812 

0325 

0369 

28.8 

7.6 

16 

0.78 

0.783 

0312 

22 

10 

AN 

sit slot 

0.689 

0.422 

0386 

0308 

0345 

0.37 

32.9 

5.8 

16 

0.7 

0382 

0.462 

23 

12 

AN 

sit slot 

0.756 

0.4 

0.434 

0334 

0374 

0377 

303 

3.7 

16 

0.68 

0.963 

0.420 

24 

14 

AN 

sit slot 

0.742 

0.41 

0.443 

0353 

0.407 

0394 

28 

33 

16 

037 

1.037 

0388 | 

30 

16 

AN 

sit slot 

0.743 

0.4 

0.466 

0366 

0.422 

0.406 

28.7 

43 

16 

0.45 

1.114 

0364 

31 

18 

AN 

sit slot 

0.772 

0J65 

0.477 

0.842 

0.427 

0.41 

313 

43 

16 

039 

1.179 

0343 

32 

20 

AN 

sit slot 

0.77 

0.356 

0.483 

0339 

0.431 

0.413 

29.6 

3.9 

16 

033 

1.240 

0325 

34 

14 

AN- Stag 

round 

0.73 

0.339 

0.422 

0.761 

0.343 

0391 

26.6 

3.6 

16 

n/a 

0.616 

0.556 

41 

22 

AN-Stag 

round 

0.779 

0383 

0386 

0.769 

0327 

0373 

41.6 

53 

16 

n/a 

0.776 

0.446 

42 

30 

AN-Stag 

round 

0.709 

0.339 

0.433 

0.772 

0353 

0352 

14 _ il 

39.8 

l 

93 

16 

n/a 

0.904 

0381 


Table-5. Mhrfng uniformity and NOx production results for the annular mixers 
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HO 

#eq 

Mixer 

Hole 

Cum Vol 

B(+) 

B(-) 

SUM-B 

AMIX 

MMIX 

Cumul 

El NOx 

Do 

Pen r/R 

Blockage 

T.E. x/R 


holes 

Type 

type 

@ f stoi 






NOx, ppm 

Global 

inch 

from Do 

mSmm 

T.E. 2x/H 

44 

6 

cir 

alg slot 

0.876 

0.199 

0.425 

0.624 

0351 

0371 

45.1 

5.97 

3 

IkEHI] 

wmM 

0.776 

51 

10 

cir 

alg slot 

0.89 | 

■ill m 

0.406 

0.596 

0357 

0356 

44.5 

KSail 

3 

0.82 

0239 

0.600 

52 

14 

cm 

alg slot 

0.87 


0364 

0.545 

0.333 

0.335 

65 

10.7 

3 

0.67 

0202 

0.507 

53 

18 

cm 

alg slot 

0.803 

0277 

0351 

0.628 

0324 

0344 

49.7 

83 

3 

0.6 

0321 

0.448 

54 

22 

cm 

alg slot 

0.794 

0.41 

0364 

0.774 

0.334 

0367 

49.1 

8.07 

3 

0.54 

0355 

0.405 

46 

26 

cm 

alg slot 

0.797 

0.409 

0377 

0.786 

0345 

0381 

37.7 

7.1 

3 

0.49 

0.386 

0.373 

55 

30 

CIR 

alg slot 

0.788 

0.423 

0392 

0.815 

0364 

0399 

38.5 

4.76 

3 

0.46 

0.414 

0.347 

56 

34 

cm 

alg slot 

0.791 

0.432 

0.403 

0335 

0.372 

0.404 

35.8 

633 

3 

0.44 

0.440 

0.325 

57 

38 

CIR 

alg slot 

0.795 

0.413 

0.404 

0.817 

0375 

0.41 

40.1 

10.1 

3 

0.44 

0.468 

0.309 

6 

6 

CIR 

round 

0.8603 

0.133 

0395 

0.528 

0.336 

0358 

42.9 

8 

3 

0.85 

0.407 

0.426 

7 

8 

CIR 

round 

0.847 

0.154 

0386 

0.54 

0.331 

0.36 

53.6 

11.4 

3 

0.72 

0.470 

0.369 

8 

10 

cm 

round 

0.8405 

0.172 

0369 

0341 

0319 

0362 

51.7 

12.6 

3 

0.66 

0.525 

0.330 

9 

12 

CIR 

round 

0.852 

0.309 

0351 

0.66 

0315 

0351 

57 

83 

3 

0.6 

0.575 

0.301 

10 

14 

CIR 

round 

0.861 

0.383 

0.342 

0.725 

0319 

0355 

52.2 

5.9 

3 

0.594 

0.621 

0.279 

16 

16 

CIR 

round 

0.852 

0.37 

0.347 

0.717 

0318 

0362 

37.9 

4.7 

3 

0.556 

0.664 

0.261 

17 

18 

CIR 

round 

0.8381 

0.392 

0356 

0.748 

0341 

0.37 

26.9 

3.9 

3 

0.543 

0.704 

0.246 

18 

20 

CIR 

round 

0.826 

0.41 

0.367 

0.777 

0.356 

0.381 

27 

23 

3 

0.523 

0.743 

0.233 

38 

24 

CIR 

round 

0.816 

0.387 

0399 

0.786 

0383 

0.39 

25.1 

7.9 

3 

0.482 

0.813 

0.213 

1 

6 

CIR 

sit slot 

0.742 

0.326 

0364 

0.69 

0321 

0.368 

37.1 

53 

3 

0.69 

0.685 

0.297 

2 

8 

CIR 

sit slot 

0.796 

0397 

0333 

0.73 

0313 

033 

48.2 

4.2 

3 

0.402 

0.790 

0.257 

3 

10 

CIR 

sit slot 

0.811 

0.4 

0375 

0.775 

0372 

0368 

25.7 

33 

3 

0.386 

0.884 

0.230 

4 

12 

CIR 

sit slot 

0.816 

0.388 

0.462 

0.885 

0.414 

0.401 

28.4 

4.6 

3 

0.358 

0.968 

0.210 

5 

14 

CIR 

sit slot 

0.813 

0.383 

0.487 

0.87 

0.431 

0.413 

32.7 

63 

3 

0.31 

1.046 

0.194 

35 

16 

CIR 

sit slot 

0.824 

0.378 

0.512 

0.89 

0.445 

0.417 

30.7 

73 

3 

0.289 

1.118 

0.182 

36 

18 

CIR 

sit slot 

0.821 

0.37 

0.514 

0.884 

0.446 

0.405 

32.9 

73 

3 

0.224 

5270 

0.171 

37 

20 

CIR 

sit slot 

0.821 

0.4 

0313 

0.913 

0.451 

0.41 

34.1 

9.2 

3 

0.244 

1.251 

0.163 


Table-6. Mixing uniformity and NOx production results for the circular mixer 
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Figure- 1 . a) Geometric configuration of the mixing section, and 

b) orthogonal view of the computational grid with a 45° slot 
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Figure- 2. 


a) Geometric configuration of the typical annular section, and 

b) orthogonal view of a single sector of the computational domain with round holes 
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Figure-3. Description of the volume histogram statistical analysis for the mixing non-uniformity evaluation 



jet penetration 2r!H from the outer wall 



Figure-4. Effect of the variation of the number of orifices on the jet penetration 
for the annular mixers 
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Figure-6. Effect of the variation of the number of orifices on the mixing non-uniformity AMIX 
for the annular mixers 




AMIX (area weighted mixing non-uniformity deviation) 





MMIX (area weighted mixing non-uniformity deviation) 



Figure-8. Effect of the variation of the number of orifices on the mixing non-uniformity MMIX 
for the annular mixers 





MMIX (area weighted mixing non-uniformity deviation) 



Figure-9. Effect of the variation of the number of orifices on the mixing non-uniformity MMIX 
for the circular mixer 
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Figure-10a. Correlation between AMIX and jet penetration far the annular mixers 
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Figure-10b. Correlation between AMIX and jet penetration far the circular mixers 


NAS A/CR— 2000-2 1 0672 









MMIX MMIX 


0.50 r 


0.45 


- 


19 ANN- AS 


▲ 

19ANN-RO 

" 

▼ 

16ANN-RO 

- 

— o— 

16ANN-SS 

- 




y~ — □-- - 


0.40 


0.35 


0.30 



19ANN =19 in D, annular mixer configuration 
16 ANN -16in D c annular mixer configuration 


AS = flow aligned slot (0 ), UW=4 

RO = round orifice 

SS = 67 . 5 “ sltanted slot UW=4 


X 


I I „.J_ 


0.4 


0.5 0.6 0.7 0.8 

jet penetration 2r!H from, the outer wall 


0.9 


1.0 


Figure-lla. Correlation between MMIX and jet penetration for the annular mixers 



Figure-lib. Correlation between MMIX and jet penetration for the circular mixers 
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19 ANN =19in D 0 annular mixer configuration 
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Figure-12. Effect of the variation of the number of orifices on themixing non-uniformity volume histogram shape 
for the annular mixers 
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Figure-13. Effect of the variation of the number of orifices on themixing non-uniformity volume histogram shape 
for the circular mixers 





Figure-14b. Correlation between SUM-B and jet penetration for the circular mixer 



16ANN-RO 









NAS A/CR— 2000-2 1 0672 


80 



Number of equivalent orifices per row 


Figure-15. Effect of the variation of the number of orifices on the mixer NO x production 
for the annular mixers 





NAS A/CR— 2000-2 1 0672 



Figure-16. Effect of the variation of the number of orifices on the mixer NOx production 
for the circular mixers 
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Figure-17a. Correlation between CUM-NOx and jet penetration for the annular mixers 



Figure-17b. Correlation between CUM NOx and jet penetration for the circular mixer 
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Figure-18. Effect of the variation of the number of orifices on the global NOx production 
for the annular mixers 
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Figure-19. Effect of the variation of the number of orifices on the global NO x production 
for the circular mixers 








Figure-20. Mixing modes for the circular (a) and annular (b) mixer types 
with representative round orifices 
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Figure-21a. Comparison between the performance of the annular and circular mixers with the aligned slot 



Figure-2 lb. Comparison between the performance of the annular mixer with round and slanted slot orifices 
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Figure-22a. Correlation between the NOx production and mixing non-uniformity (SUM-B) 
for the annular mixers (phase-b) trend 
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Figure-22b. Correlation between the NOx production and mixing non-uniformity (SUM-B) 
for the circular mixers (phase-b) trend 


NAS A/CR— 2000-2 1 0672 







BALEINOx 



Figure-2 3a. Correlation between the NOx production and mixing non-uniformity (AMIX) 
for the annular mixers (phase-b) trend 
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Figure-23b. Correlation between the NOx production and mixing non-uniformity (AMIX) 
for the circular mixers (phase-b) trend 
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Figure-24a. Correlation between the NOx production and mixing non-uniformity (MMIX) 
for the annular mixers (phase-b) trend 
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MMIX 

Figure-24b. Correlation between the NOx production and mixing non-uniformity (MMIX) 
for the circular mixers (phase-b) trend 
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Figure- Al. Configuration HO- 1 , 6 7. 5° slot UW= 4, 
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Figure- A2. Configuration HO -2, 67.5 0 slot L/W— 4, 8 orifices /row 
3=52.0, MR=2.96, DR=2.28, ty RZ = 1.80, = 0.416 ( 2-adj acent sectors) 
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Figure- A3. Configuration HO-3 , 67.5° slot LTW-4, 
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Figure- A4. Configuration 110-4, 67.5° slot LfW=4, 12 orifices /row 
J=52.0, MR- 2.96, DR=2.2S , § RZ = 1.80, ^ = 0.416 ( 2-adjacent sectors) 
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Figure-AS. Configuration HO~5, 67.5° slot L/W=4, 14 orifices /row 
J=52.0, MR=2.96, DR=2.28 , (j)^ = 1.80, ^ = 0.416 ( 2-adj acent sectors) 
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Figure- A6. Configuration HO-6 , 6 round orifices /row 


J=52.0, MR-2.96, DR-2.28, 0 fiZ = 1.80, = 0.416 ( 2-adj acent sectors) 
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Figure- A7. Configuration HO-7, 8 round orifices / row 


J=S2.0, MR= 2.96, DR- 2.28, § RZ a 1.80, = 0.416 ( 2-adj acent sectors) 
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Figure-A8 . Configuration HO -8, 10 round orifices /row 


J=52.0 , MR-2.96, DR=2.28, § KZ - 1.80, § lz - 0.416 ( 2-adj acent sectors) 
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Figure -A9. Configuration HO-9, 12 round orifices /row 


J-52.0, MR- 2.9 6, DR- 2.28, § RZ = 1.80, = 0.416 ( 2-adjacent sectors ) 
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Figure- A10. Configuration HO- 10, 14 round orifices /row 


J-52.0, MR=2.96, DR=2.28, § RZ = 1.80 , ^ LZ = 0.416 ( 2-adjacent sectors) 
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Figure-All. Annular mixer configuration HO- 11, 29 round holes per row (6 eq holes/row) 
J=52.0, MR=2.96 , DR=2.28, 1.80, 0.416 (2 sectors) 
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Figure-Al2. Annular mixer configuration HO -12, 39 round holes per row (8 eq holes /row) 
J-52,0 f MR=2.96, DR=2.28, 4 ^= 1.80 , 4 ^= 0.4 16 (2 sectors) 
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Figure-A13. Annular mixer configuration 110-13, 48 round holes per row (10 eq holes/row) 
J=52.Q, MR— 2.96, DR=2.28 t ^1.80, ^ 0.416 (2 sectors) 
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Figure-A14. Annular mixer configuration HO-14 , 58 round holes per row (12 eq holes/row) 
J-52.0, MR=2.96, DR=2.28, ^ 1.80 , §^=0.416 (2 sectors) 
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Figure-AlS . Annular mixer configuration HO- 15, 68 round holes per row (14 eq holes /row) 
J=52.0, MR=2.96 t DR=2.28, ^=1.80, ^0.416 (2 sectors) 
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Configuration 110-16, 16 round orifices /row 


1—52.0, MR— 2,96, DR— 2.28, - 1.80, = 0.416 ( 2-adj acent sectors) 
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Figure- Al 7 . 
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Configuration HO-1 7, 18 round orifices /row 


J-52.0, MR=2,96, DR=2.28, § RZ = 1.80, - 0.416 ( 2-adjacent sectors) 
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Figure- AI8. Configuration HO- 18, 20 round orifices /row 


J—52.0 } MR- 2.96, DR-2.28, § R7 = 1.80, - 0.416 ( 2-adj aceni sectors) 
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Figure- A19. Annular mixer configuration HO-1 9, 19 round holes per row (4 eq holes fro w ) 
J-52.0, MR= 2.96, DR= 2.28, ^ 1.80, ^ 0.416 (2 sectors) 
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Figure-A20. Annular mixer configuration HO-20 , 67.5° slanted slots UW—4 , 29 per row (6 eq holes/ row) 
J=52.0, MR- 2.96, DR=2.28, ^ 1.80 , ^0.416 (2 sectors) 
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Figure-A21 , Annular mixer configuration HO-21 , 67.5 9 slanted slots LfW= 4, 39 per row (8 eq holes/ row) 
J-52.0, MR=2.96, DR=2.28, §^=1.80, <j ^-0.416 (2 sectors) 
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Figure- A22. Annular mixer configuration HO -22, 67.5 0 slanted slots LTW=4, 48 per row (10 eq holes/row) 
J=52.0, MR=2.96, DR— 2.28, §^=1.80, §^=0.416 (2 sectors) 
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Figure-A23. Annular mixer configuration BO-23, 67.5° slanted slots L/W-4, 58 per row (12 eq holes /row) 
J-52.0, MR=2.96, DR-2.28, ^ 1.80 , 0.416 (2 sectors) 
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Figure-A24. Annular mixer configuration HO-24, 67.5 0 slanted slots UW- 4, 68 per row (14 eq holes/row) 
J=52.0, MR=2.96, DR-2.28 , ^=1.80, ^=0.416 (2 sectors) 
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Figure- A2 5. Annular mixer configuration HO-25, 77 round holes per row (16 eq holes/row) 
7=52.0, MR=2.96 t DR-2.28, §^=1.80, §^=0.416 (2 sectors) 
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Figure -A2 6. Annular mixer configuration HO-26, 87 round holes per row (18 eq holes/row) 
J-52.0, MR=2.96, DR=2.28, ^1.80, ^=0.416 (2 sectors) 
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Figure-A27. Annular mixer configuration HO-27, 97 round holes per row (20 eq holes/row) 
J=52.0, MR=2.96 , DR=2.28, §^=1.80, ^=0.416 (2 sectors) 
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Figure -A2 8. Annular mixer configuration HO -28, 106 round holes per row (22 eq holes/row) 
7=52.0, MR=2.96 , DR=2.28, §^=1.80, § u =0.416 (2 sectors) 


NAS A/CR— 2000-2 1 0672 


00 

o 



x/H- 0,5 


x/H= 0.37S\ 


j /H=0.2S 


TEMP , K 

2500 
2300 
2100 
1900 
1700 
1500 
1300 
1100 
900 




.«///= 0.5 


x/H= 0.375 


x/H=0.25\ 

NO, 

log(gr NQJ {sec m 3 J) 




x/H =0.5 


x M-Q.37S 


x/II=0.25 


Normalized (j) 


x/H=0.125 


x/H=0.08 



Figure-A29. Annular mixer configuration HO-29 , 116 rouiui holes per row (24 eg holes/row) 
J=S2.0, MR-2.96, DR-2.28, §^=1.80, § lJZ =0.416 (2 sectors) 
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Figure-A30 ♦ Annular mixer configuration HO -30, 67.5° slanted slot , L/W=4, 77 per row (16 eq holes/row) 
J-52.0, MR- 2. 96, DR=2.28, 1.80 , $^=0,416 (2 sectors) 
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Figure-A31. Annular mixer configuration HO-31, 67. 5° slanted slot, L/W-4 87 per row (18 eq holes/row) 
7-52.0, MR=2.96, DR=2.28, fes 1.80 , ^=0.416 (2 sectors) 
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Figure-A32. Annular mixer configuration HO~32 s 67.5 0 slanted slot , L/W~4 f 97 per row (20 eq holes/row) 
J=52.0, MR- 2.96, DR=2.28, 0 ^= 1.80 , ^-0.416 (2 sectors) 
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Figure-A33. Annular mixer configuration HO-33. 145 round holes per row (30 eq holes/row) 
7=52.0, MR=2.96, DR=2.28, ^ 1.80 , 0.41 6 (2 sectors) 





NAS A/CR— 2000-2 1 0672 


00 

Lh 



0.125 


TEMP , K 

2500 
2300 
2100 
1900 
1700 
1500 
1300 
1100 
900 




Normalized <|> 


mam 


Figure-A34. Annular mixer configuration HO-34, 68 round holes per row (14 eq holes/row) 
J=52.0, MR-2.96, DR=2.28, ^ 1.80 , ^ 0.416 (2 sectors) [STAGGERED HOLE CONF] 
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Figure- A3 5. Configuration HO-35 , 67.5 0 slot L/W— 4, 16 orifices /row 


J=52.0 , MR- 2.96, DR- 2.28, = 1.80, 0 ^= 0.416 ( 4-adjacent sectors) 
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Figure- A3 6. 
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Configuration HO-36, 67.5° slot lfW-4, 18 orifices /row 


J=52.0, MR-2,96, DR=2.28 , § RZ - 1.80, (J)^ = 0.416 ( 4-adjacent sectors) 
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Figure- A3 8. Configuration HO-38, 24 round orifices l row 


J-52.0, MR- 2.9 6, DR=2,28, § RZ = 1,80, ^ = 0.416 ( 6-adjacent sectors) 
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Figure -A3 9. Annular mixer configuration HO-39, 64 round holes per row (6 eq holes/ row) 
J=52.0, MR— 2.96, DR=2.28, <j ) R2 —1.8Q, 0.416 (2 sectors) 
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Figure-A40. Annular mixer configuration HO-40, 21 round holes per row (2 eq holes f row) 
J=52.0 f MR=2.96, DR=2.28, ^=1.80, 0.416 (2 sectors) 
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Figure-A41. Annular mixer configuration HO-41 , 106 round holes per row (22 eq holes/row ) 
J=S2.0, MR- 2.96, DR=2.28 , 1.80 , <^= 0.416 (2 sectors) [STAGGERED HOLE CONE] 
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Figure- A42, Annular mixer configuration HO-42, 145 round holes per row (30 eq holes/row) 
J=52.0, MR=2.96, DR=2.28, $ K = 1.80 , $^=0,416 (2 sectors) [STAGGERED HOLE CO NF] 
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Figure-A43. Annular mixer configuration HO-43,flow aligned slot LSW-4, 64 per row (6 eq holes/row) 
J=52.0, MR=2.96, DR=2.28, fas 1.80 , ^=0.416 (2 sectors) 
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Figure-A44. Configuration HO -44, flow aligned slot LPN-4 , 6 orifices /row 
J-S2.0, MR- 2.9 6, DR-2,28 t 0^ = 1.80, 0 /j; = 0.416 ( 2-adjacent sectors) 
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Figure-A45. Annular mixer configuration HO -45, flow aligned slot LTW-4, 277 per row (26 eq holes/row) 
J=52.0, MR=2.96, DR=2.28, §^=1.80, §^=0.4 16 (10 sectors) 
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Configuration HO-46, flow aligned slot L/W=4, 26 orifices /row 
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Figure- A47 . Annular mixer configuration HO -47, flow aligned slot UVf—4, 167 per row (10 eq holes/row) 
J=52.0, MR- 2.96, DR=2.28, <j >^=1.80, §^=0.416 (5 sectors) 
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Figure-A48. Annular mixer configuration ftO~48, flow aligned slot L/W=4, 149 per row (14 eq holes/row) 
J=52.0 , MR=2.96 > DR=2,28 } ^ =1.80, ^=0.416 (5 sectors) 
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Figure -A49. Annular mixer configuration HO -49, flow aligned slot UW=4, 192 per row (18 eq holes/row) 
J=52.0, MR- 2.96, DR=2.28, ^ 1.80 , ^ 0.416 (6 sectors) 
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Figure- A50. Annular mixer configuration HO -50, flow aligned slot L/W=4 f 235 per row (22 eq holes/ row) 
J=52.0, MR-2.96, DR-2.28, §^=1.80, §^=0.416 (5 sectors) 
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Figure- A52. Configuration HO-52, flow aligned slot UW- 4, 14 orifices /row 


J=52.0 MR-2.96, DR=2.28, § KZ = 1.80, - 0.416 ( 4-adjacent sectors) 
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Figure~A53. Configuration HO-53, flow aligned slot LfW-4, 18 orifices /row 


J-52.0, MR-2.96, DR- 2.28, (j)^ = 1,80, ^ LZ = 0.416 { 4-adjacent sectors) 
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Figure-A54. Configuration HO-54, flow aligned slot UW-4, 22 orifices l row 


J=52.0 , MR- 2.96, DR= 2.28, <\> RZ - 1.80, 0^ = 0.41 6 ( 6-adjacent sectors ) 
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Figure-ASS. Configuration HO-55, flow aligned slot L/W-4, 30 orifices /row 


J-S2.0 , MR- 2.9 6, DR- 2. 28, § RZ = 1.80, - 0.416 ( 7 •adjacent sectors) 
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Fig ure -AS 6. Configuration HO- 5 6, flow aligned slot UW— 4, 34 orifices l row 
J=52.0, MR— 2.96, DR=2.28, § RZ — 1.80, = 0.416 ( 8 -adjacent sectors) 
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Figure- A5 7. Configuration HO-57, flow aligned slot L/W=4, 38 orifices /row 
J-52.0 , MR- 2.96, DR-2.28, 4 ^ = 1.80, 0 ^ = 0.416 ( 9 -adjacent sectors) 
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Figure-A58. Annular mixer configuration HO-58, flow aligned slot LfW=4, 299 per row (28 eq holes/row) 
J-52.0, MR=2.96, DR=2.28, ^=1.80, ^=0.416 (10 sectors) 
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Figure -A59. Annular mixer configuration HO-59, flow aligned slot UW^4, 21 per row (2 eq holes /row) 
J-52.0, MR- 2.96, DR=2.28, ^=1.80, §^=0.416 (1 sector) 
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Figure -A60. Annular mixer configuration HO -60, flow aligned slot UW-4 , 11 per row (1 eq holes/row) 

MR=2.96 t DR=2.28, $^=1.80, ^=0,416 (I sector) 
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b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-Bl. Configuration HO-1. Circular mixer, 6 round holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B2. Configuration HO-2. Circular mixer, 67.5° slot L/W=4, 8 holes/row 
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a. Planar deviations (AMIX and MMIX ) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-83. Configuration HO -3. Circular mixer, 67.5° slot UW= 4, 10 holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B4. Configuration HO -4. Circular mixer, 67.5° slot L/W=4, 12 holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B5. Configuration HO-5. Circular mixer, 67.5° slot L/W=4, 14 holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B6. Configuration HO-6. Circular mixer, 6 round holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B7. Configuration HO -7. Circular mixer, 8 round holes/row 
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a. Planar deviations (AMIX and MMIX ) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 


Figure-B8. Configuration HO-8. Circular mixer , 10 round holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 


Figure-B9. Configuration HO-9. Circular mixer , 12 round holes/row 
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a. Planar deviations ( AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B10. Configuration HO -10. Circular mixer, 14 round holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-Bll. Configuration HO-11. Annular mixer, 6 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B12. Configuration HO-12. Annular mixer, 8 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b . mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B13. Configuration HO-13. Annular mixer , 10 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x!H < 1 


Figure-B14. Configuration HO-14. Annular mixer, 12 round equivalent holes/row 
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a . Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B15. Configuration HO -15. Annular mixer, 14 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B16. Configuration HO -16. Circular mixer, 16 round holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B17. Configuration HO -17. Circular mixer, 18 round holes/row 


NAS A/CR— 2000-2 1 0672 


128 


Normalized volume fraction per bin 





Normalized cumulative volume fraction 




a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B18. Configuration HO-18. Circular mixer, 20 round holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B19. Configuration HO -19. Annular mixer, 4 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX ) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 


Figure-B20. Configuration HO-20. Annular mixer , 67.5° slot L/W=4, 6 equiv holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 


Figure-B21. Configuration HO-21. Annular mixer, 67.5° slot UW-4, 8 equiv holes/row 
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a . Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B22. Configuration HO-22. Annular mixer, 67.5 0 slot UW=4, 10 equiv holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B23. Configuration HO-23. Annular mixer, 67.5° slot L/W=4, 12 equiv holes/row 
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a . Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 


Figure-B24, Configuration HO-24. Annular mixer > 67.5° slot L/W—4, 14 equiv holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B25. Configuration HO -25. Annular mixer, 16 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-1126. Configuration HO-26. Annular mixer > 18 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B27. Configuration HO-27. Annular mixer, 20 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 


Figure-B28. Configuration HO-28. Annular mixer , 22 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX ) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 


Figure-B29. Configuration HO-29. Annular mixer , 24 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX ) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B30. Configuration HO-30. Annular mixer, 67.5° slot L/W=4, 16 equiv holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B31. Configuration HO-31. Annular mixer , 67.5° slot UW=4, 18 equiv holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B32. Configuration HO-32. Annular mixer, 67.5 0 slot UW—4, 20 equiv holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B33. Configuration HO-33. Annular mixer, 30 round equivalent holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B34. Configuration HO -3 4. Ann-staggered mixer, 14 round equiv holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B35. Configuration HO-35. Circular mixer, 67.5° slot UW=4, 16 holes/row 
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a. Planar deviations ( AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B36. Configuration HO -3 6. Circular mixer, 67.5° slot UW= 4, 18 holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B37. Configuration HO-37. Circular mixer, 67.5 0 slot UVf-4, 20 holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B38. Configuration HO-38. Circular mixer , 24 round holes/row 
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a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B39. Configuration HO-39. Annular mixer, 6 round equivalent holes/row 
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Axial distance 2x/H from the leading edge of the orifice 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 


Figure-B40. Configuration HO-40. Annular mixer, 2 round equivalent holes/row 
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Normalized volume fraction per bin 






Axial distance 2x/H from the leading edge of the orifice 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B41. Configuration HO-41. Ann-staggered mixer, 22 round equivalent holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 



b. 
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MAINSTREAM 


Non-dimensional equivalence ratio (f) 

a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B42. Configuration HO-42. Ann-staggered mixer , 30 round equivalent holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 



Axial distance 2x/H from the leading edge of the orifice 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B43 . Configuration HO-43. Annular mixer Jlow aligned slot L/W=4, 6 eq holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 


a. 



Anal distance x/R from the leading edge of the orifice 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 


Figure-B44. Configuration HO-44. Circular mixer, flow aligned slot L/W=4, 6 holes/row 
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Normalized volume fraction per bin 





a. 



Axial distance 2x/H from the leading edge of the orifice 


b. 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 


Figure-B45. Configuration HO-45. Annular mixer Jlow aligned slot L/W=4, 26 eq holes/row 
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Normalized volume fraction per bin 







a. Planar deviations (AMIX and MMIX ) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B46. Configuration HO-46. Circular mixer, flow aligned slot L/W=4, 26 holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 



Axial distance 2x/H from the leading edge of the orifice 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B47. Configuration HO-47. Annular mixer Jlow aligned slot L/W=4, 10 eq holes/row 
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Normalized volume fraction per bin 






b. 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B48. Configuration HO -48. Annular mixer Jlow aligned slot UW=4, 14 eq holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 


a 



Axial distance 2x/H from the leading edge of the orifice 



a. Planar deviations (. AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B49. Configuration HO -49. Annular mixer Jlow aligned slot L/W=4, 18 eq holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 0 Mixing non-uniformity 



Axial distance 2x/H from the leading edge of the orifice 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figwe-B50. Configuration HO-50. Annular mixer Jlow aligned slot UW=4, 22 eq holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 



b. 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B51. Configuration HO-51. Circular mixer, flow aligned slot L/W=4, 10 holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 


a. 



Axial distance xJR from the leading edge of the orifice 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 


Figure-B52. Configuration HO-52. Circular mixer, flow aligned slot L/W=4, 14 holes/row 
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Normalized volume fraction per bin 






Axial distance x/R from the leading edge of the orifice 


b. 



a. Planar deviations (AMIX and MM1X ) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 


Figure-353. Configuration HO -53. Circular mixer, flow aligned slot UW-4, 18 holes/row 


NAS A/CR— 2000-2 1 0672 


164 


Normalized volume fraction per bin 





Normalized cumulative volume fraction 


a. 



Axial distance xJR front the leading edge of the orifice 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 


Figure-B54. Configuration HO -54. Circular mixer, flow aligned slot L/W=4, 22 holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 




a. Planar deviations (AMIXand MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B55. Configuration HO-55. Circular mixer, flow aligned slot LfW=4, 30 holes/row 
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Normalized volume fraction per bin 





a. 



Axial distance x/R from the leading edge of the orifice 


b. 



a. Planar deviations (AMIX and MMIX ) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 

Figure-B56. Configuration HO-56. Circular mixer, flow aligned slot L/W=4, 34 holes/row 
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Normalized volume fraction per bin 





Normalized cumulative volume fraction 



b. 



Non-dimensional equivalence ratio (f) 

a. Planar deviations ( AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < x/R < 1 
Figure-B57. Configuration HO-57. Circular mixer, flow aligned slot L/W=4, 38 holes/row 
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Normalized volume fraction per bin 
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Axial distance 2x/H from the leading edge of the orifice 

b. 



a. Planar deviations (AMIX and MMIX ) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B58. Configuration HO -5 8. Annular mixer Jlow aligned slot L/W=4, 28 eq holes/row 
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Normalized volume fraction per bin 






Axial distance 2x/H from the leading edge of the orifice 

b. 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B59. Configuration HO-59. Annular mixer Jlow aligned slot L/W=4, 2 eq holes/row 
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Normalized volume fraction per bin 





Axial distance 2x/H from the leading edge of the orifice 



a. Planar deviations (AMIX and MMIX) results 

b. mixing non-uniformity volume histograms for 0 < 2x/H < 1 

Figure-B60. Configuration HO -60. Annular mixer flow aligned slot UW= 4, 1 eq holes/row 
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Normalized volume fraction per bin 
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